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ABSTRACT 

We present near-IR JH spectra of the central regions of the dwarf starburst galaxy 
NGC 1569 using the Florida Image Sheer for Infrared Cosmology and Astrophysics 
(FISICA) . The dust-penetrating properties and available spectral features of the near- 
IR, combined with the integral field unit (IFU) capability to take spectra of a field, 
make FISICA an ideal tool for this work. We use the prominent [Hei] (1.083/im) and 
Pa/3 (1.282 (im) lines to probe the dense star forming regions as well as characterize 
the general star forming environment around the super star clusters (SSCs) in NGC 
1569. We find [Hei] coincident with CO clouds to the north and west of the SSCs, 
which provides the first, conclusive evidence for embedded star clusters here. 

Key words: galaxies: active — galaxies: dwarf — galaxies: star clusters infrared: 
galaxies 



1 INTRODUCTION 



NGC 1569 is a nearby, 2.2 Mpc (|lsrael fc van Driell Il99cl ) . 
dwarf starburst galaxy. Its centre contains two, prominent 
super star clusters (SSCs) , designated "A" and "B" by 
Ide Vaucouleurs et alJ (ll974Tl . Originally thought to be fore- 
ground stars ( Able j|l971 ) . later spectroscopic studies con- 
cluded th at they are s tar clusters a s sociat ed with NGC 
1569 (cf. |Pradafll994h . IHunter et ail (|2000h performed a 
high- resolution, HST optical study of the stellar populations 
around the SSCs. They fit ages of ^7 Myr for A and 10-20 
Myr for B. 

This galaxy has undergone star formation for almost 
a complete Hubble ti me, containing clusters spanning ages 
from 4 Myr to 10 Gyr lUloisi et al. l l200ll : lAnders et alj|2004 
IHunter et a"i1l2000l ). IWallerl (|l99lT) sug gests this syst e m un - 
derwent six starbursts in its history. IHunter et alJ (|2000l ) 
find most clusters are <30 Myr, indicating most cluster for- 
mation occurr e d at the end of the last starburst episode. 
iPasquali et alJ (1201 ll ) point out that each consecutive star- 
burst has a shorter duration and a higher star formation 
rate. Radio observations show a ~200 p c diameter hole in 
the gas distribution ar ound the SSC A (|lsrael fc van Driel 
ll990l : lAloisi et al.ll200ll ). This indicates a lack of star forma- 



tion in this region around the SSCs. This inne r region is also 
highly disturbed l|Westmoquette et al.l 200711) and lacks su- 
pernova remnants (SNRs) l|Greve et al 2002). Furthermore, 



the more massive clusters formed at earlier ti mes, while at 
later t imes, only less mass ive clusters formed (|Anders et al.l 
l2004h . I Anders et all l|2004l) attribute this to the strong, clus- 
ter winds of the massive clusters, which prevented large 
molecular clouds from forming. 

Star f ormation may still be occurring in a region t o 
the north ijWestmoauette et al.ll2q07bl: iTokura et al]l2006h . 
near cluster 30 as designated by IHunter et al.l (l2000h an d 
to the north west, around cluster 10 ( Hunter et al. 2000f ). 
The north west region is quite re markable. It is domi- 
nated by cluster 10, designation from IHunter et al", I (I2000T) 
which was initially thought to be another SSC ( Pradal 
Il994h. but later evidence suggests it is a less massi ve cluster 
i Westmoauette et al.ll2007al : iBuckalew et al.ll2000l ). This re- 
gion itself contains many clusters, with ages of ~25 Myr 
l| Anders et alJ 120041). Just south of 10 lies an Hn region 
l|Buckalew et al.l I2000T) and to it's west extends a giant 
molecul ar cloud (I Taylor et al.l Il999t ) . Using mid-IR obser- 
vations, ITokura et alJ (|2006l ) find possible evidence for an 
embedded cluster here, which could substanti ate evidence 
that star formation is progressing to the west (|Greve et al.1 
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2002). |Pasquali et alj (|201 lh find that the dust extinction is 
higher in the NW than in the SE. 

In this article we use the Florida Image Slicer for In- 
frared Cosmology and Astrophysics (FISICA), a near-IR 
image-slicing integral field unit (IFU), to observe the en- 
vironment around the SSCs and cluster 10. Near-IR wave- 
lengths are particularly useful as they contain the [He i] line 
at 1.083 fim and the Pa/3 line at 1.282 /im, both of which 
are diagnostic spectral features of massive stars and star for- 
mation. In addition, the dust-penetrating properties of IR 
wavelengths allows us to probe dense regions hidden at other 
wavelengths. 

Our goal in this work is to search for more evidence of 
massive star formation in the embedded regions around the 
SSCs as well as complement previous work by adding to the 
understanding of this complex environment. We present our 
work in the following manner: §2 outlines our observation 
method and analysis of the data, §3 presents a discussion of 
our results, comparing in detail to past work on this galaxy, 
and we summarize our findings in §4. 



2 OBSERVATIONS AND DATA ANALYSIS 

On October 20, 2004, we acquired JH spectra of NGC 1569 
using FISICA on the KPNO 4-m telescope. We took five 
consecutive 300 s exposures, producing a total on-source ex- 
posure of 25 min. Due to the extended nature of our target, 
we moved off-source by 95 arcsec to the SW to measure the 
sky background. 

FISICA (jEikenberrv et all I2006T ) has a field-of-view 
(FOV) of 16 by 33 arcsec at f/15. The IFU divides the field 
into 22 strips ~0.7 arcsec wide by 33 arcsec long. These 
strips are rearranged in a line and sent as a pseudo-longslit 
through the Florida Multi-object Imaging Near-IR Grism 
Observational Spectrometer (FLAMINGOS) spectrograph. 
All 22 spectra are imaged on the FLAMINGOS, 2Kx2K 
HgCdTe detector, and each has a resolution of R ~1300. 

We reduced these observations using the data pipeline, 
Florida Analysis Tool Born Of Yearning for high quality 
scientific data (fatboy). fatboy (C. Warner, et al. 2012, in 
preparation) calibrates each spectral frame using flat fields 
and dark frames and corrects for cosmic rays, and bad pix- 
els. All spectra are rectified and all spectral lines are aligned 
using skylines. We did not need to worry about image distor- 
tion as a function of wavelength, because spatial information 
is preserved for IFUs that use mirrors. 

A useful tool for analyzing IFU data is a data cube. To 
make the data cube, we used a self-written Python program. 
The main difference between a standard image is each pixel 
has spectral information. By this information in the z direc- 
tion, we produced a data cube for the FISICA field of NGC 
1569. To ensure Nyquist sampling, the image of each slice 
is spread across two pixels. Considering we had a seeing of 
~3.0 pixels, we are Nyquist sampling the Airy disk and this 
is appropriate. Therefore, the reproduced image of the field 
will be foreshortened by a half in the x direction. Artificially 
spreading the flux from each slice across a width of two pix- 
els, we produced an image with the correct aspect ratio, 
having pixels with dimensions of ~0.3 arcsec. More specifi- 
cally, we measured the centroid in flux along each slice. From 
this, we measured the fractional flux for either half of the 




Figure 1. False colour image of the central regions of NGC 1569. 
North is to the upper-left and east is to the lower-left. Green 
is continuum emission, red is Pa/3 emission and blue is [Hel] 
emission. The Pa/3 and [Hel] emission almost coincide spatially. 
Diamonds designate Hell emission peaks, circles [Om] emission 
peaks, crosses are [Siv] emission peaks, and squares are positions 
of H2 regions (see text for details). The thick, white cross is the 
peak in the [He 1] emission and the thick, white 'X' is the peak in 
the Pa/3 emission for their respective images. See §3 for a discus- 
sion of the various emission peaks. 

slice. This fractional flux was then placed in each of the two 
pixels for the newly binned slice. 

2.1 Spectral Maps 

When we initially examined the spectra, we immediately 
noticed a prominent [Hel] line in many of the slices. This 
type of emission is often associated with massive stars and 
star formation. Considering the importance of this emission, 
we decided to explore its extent and relation to the SSCs in 
NGC 1569. 

An easy way to address this query is to make a 
continuum-subtracted [Hel] map. Before making this map, 
we needed to measure the width of the [He 1] line and deter- 
mine the wavelength range that would include this emission. 
First, we summed up all spectra from all slices, producing 
a single spectrum for the entire IFU frame. Fitting a sixth 
order, Gaussian function to the [Hel] line, we measured the 
full width at half maximum (FWHM) and found a value of 
FWHM = 10.2 A. We defined an 'extraction width', w ex t 
— 12.8 A (3a of the Gaussian fit). After summing up all 
[He 1] emission within ±w ext of the [He 1] line centre for each 
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Figure 2. [Hel] equivalent width map. Darker shades correspond to regions of higher EW. Notice the large equivalent width corresponds 
with the bright emission region in the NW and the lower equivalent width around the two SSCs. The plots on the right show vertical 
and horizontal cuts through each of the four regions discussed in §2.1.1. The regions are labeled as follows: A is SSC A, B is SSC B, C 
is the region around cluster 10, and D is the [Hel] peak. 



slice, we recombined the slices, producing a [Hel] map. We 
repeated this procedure for two continuum regions on ei- 
ther side of the [Hel] line, also sampling within ±w ext of 
each position, producing a pair of continuum maps. Averag- 
ing the continuum maps, we produced a master, averaged- 
continuum map. Subtracting this from the [Hel] map, we 
made the continuum-subtracted [Hel] map. 

The continuum-subtracted [Hel] map clearly shows an 
abundance of [He i] at the top, near the set of bright clusters, 
but this emission all but disappears for the two SSCs at the 
middle and bottom of the frame (see Fig. 1). 

A second prominent line apparent in the spectrum is 
Pa/3 (1.282 fim). Like [Hel], Pa/3 is a signature of young, 
massive stars, making it a useful diagnostic for star forma- 
tion in NGC 1569. Therefore, we also produced a continuum- 
subtracted map for this line as well. Following the same 
procedure outlined above for [Hel], we stacked all spectra 
on the frame and fit a sixth order Gaussian function to the 
Pa,5 line. Measuring a FWHM of 11.3 A, we defined w ext = 
14.1 A(3<t of the Gaussian fit), extracted all Pa/3 emission 
and produced a Pa/3 map. Selecting continuum positions on 
either side of the Pa/3 line, we extracted continuum emission 
using Tiiwt, where w ex t is the same size as that used for 
the Pa/3 map. Recombining these continuum regions and av- 
eraging, we produced an averaged-continuum map. Finally, 
subtracting the averaged-continuum map from the Pa/3 map 
left us with a continuum-subtracted Pa/3 map. 



Before continuing, we performed a more quantitative 
analysis to address how much stronger, or weaker, the [Hel] 
emission is compared to the continuum for specific regions 
in the map. We chose regions with a noticeable difference in 
emission; specifically the two SSCs, A (centre; region 'A') 
and B (bottom; region 'B'), the region around cluster 10 
(region 'C') and the prominent [Hel] emission region to it's 
west, right, (region 'D') (see Fig. 2 and Fig. 3). Placing aper- 
tures over each region, we measured the total [Hel] flux in 
each aperture. We defined the aperture radius as 4.3 pixels, 
which is twice the FWHM of the Gaussian point-spread- 
function of the SSCs. This choice in aperture size was suf- 
ficient to encompass each SSC without a substantial con- 
tribution from background sources. Dividing the fluxes by 
the total flux for the entire [He i] continuum-subtracted map, 
gave us the fraction of [He i] emission each region contributes 
to the IFU frame (see Table 1). Repeating this process with 
the same-sized apertures for the averaged-continuum map 
gave us the fraction of continuum emission each region con- 
tributes to the IFU frame (see Table 1). Notice the increase 
in intensity of [Hel] flux from bottom-to-top, region 'B' to 
'A' to 'C to 'D', while the opposite effect is seen for the 
continuum flux. 

Repeating the same quantitative analysis on the Pa/3 
map that we performed on the [Hel] map, we computed 
the flux ratio in Pa/3 flux for regions A, B, C and D to 
total Pa/3 flux for the image. Lastly, using both continuum- 
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Table 1. Flux Fractions for Select Regions of NGC 1569 





Region A 


Region B 


Region C 


Region D 




(SSC A) 


(SSC B) 


(Cluster 10 Region) 


([Hel] Pk.) 


[Hel]/[Hel] toi 


3.36X10" 4 


1.03X10" 4 


1.36xl0~ 3 


1.79x 10" 3 




2.60 xicr 5 


2.21X10" 5 


3.33xl0~ 5 


6.19X10" 5 


Pa/3/Pa/3 tot 


2.27xl0~ 4 


1.31X10" 4 


1.13xl0~ 3 


1.58X10" 3 




8.98xl0~ 6 


5.44X10" 6 


5.66xl0~ 5 


1.19xl0~ 4 


Cont./Cont. tot 


2.43xl0~ 3 


3.04X10" 3 


1.21xl0~ 3 


4.82xl0~ 4 




3.22xl0~ 4 


4.84X10" 4 


4.06 xl0~ 5 


2.81xl0~ 5 


EW 


0.48 


0.13 


4.00 


18.33 




0.34 


0.18 


1.08 


7.87 



All ratios and EW are unit less. Uncertainties are listed below each 
flux ratio and EW measurement. 



subtracted maps, [Hel] and Pa/3, plus an averaged contin- 
uum map, we produced a false colour image, high-lighting 
the various emission regions (Fig. 1). 

We obtained further information on the central parts 
of NGC 1569 by making an equivalent width (EW) map. 
This is simply the [He i] continuum-subtracted map divided 
by the averaged-continuum map. Examining this EW map 
shown in Fig. 2, shows a region of large EW in the NW as 
well as a second peak to the north of cluster A. The highest 
EW is in the NW region and has a value of ~64 (See Fig. 2). 
The peak in EW to the north is about a factor of ten lower 
in EW. Also, notice the much lower EW at the positions of 
the two SSCs. Next to the EW map we show plots of vertical 
and horizontal cuts through the centre of each region. Each 
cut is the median average of ±2 pixels from the line. 



2.2 Spectra 

In addition to spectral maps, we also made spectra from 
the four regions A, B, C and D defined in §2.1 (see Fig. 3). 
These spectra were taken from the centre of each region and 
were produced by median combining the spectra ±2 pixels 
from the centre in the y-direction. Notice that the spectra 
for SSC B (Region B) is dominated by the continuum, while 
the spectra for the [Hel] peak (Region D) is dominated by 
the emission lines [Hel] and Pa/3. 



3 DISCUSSION 

Thus far we have presented a description of the [Hel] and 
Pa/3 emission in the central regions of NGC 1569. We now 
proceed by comparing our observations to other studies at 
a variety of wavelengths to try to increase understanding of 
this complex region. 

As mentioned above, there is a distinct lack of both 
[Hel] and Pa/3 emission present in both SSCs. While SSC A 
may contain trace amounts of these emission lines, SSC B 
is clearly devoid of any emission at these wavelengths. Hell 
obse rvations also show a similar trend in emission proper- 
ties (jBuckalew et al. 2000). These observations are in agree- 
ment with H I observations, which indica te a 200 pc hole in 
the H I distribution, cen t ered on SSC A (llsrael fc van Driell 
Il99d ; IStil fc Israeli 1 19981 ; lAloisi et all 1200 ll ). This suggests 
that star formation has ceased in the inner regions of NGC 
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Figure 3. Spectra through each region as labeled in Figure 2. 
Notice that the continuum emission is dominant in the spectra 
for SSC B, but the emission lines [Hel] and Pa/3 are dominant 
near the [Hel] peak. 
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1569 (|Aloisi et al.ll200lh . A detailed optical study of the stel- 
lar populations in these clusters indicates A contains lumi- 
nous blue v ariables (LBVs), w hile B is dominated by red 
supergiants (jHunter et alJl200Cf ). Thus the [Hei] emission is 
most likely from the hot, blue stars in A, while the more 
evolved, massive stars in B explains why we do not see any 
emission in this cluster. 

To the north of SSC B lies a small clump of [Hei] and 
Pa/3 emission. This region exhibits a relatively large equiva- 
lent width compared to the fie ld, peaking at ~23. It is south 
of cluster 30, as designated bv lHunter et al.l (2000), which is 
off our field. As seen in Fig. 1, there is no continuum emis- 
sion, which is consist ent with T/-band i mages, where it is 
only faintly detected dHunter et al. 2000|). as well as recent 
optical IFU observations ( Westmoquette et al. 2007b). This 
contrasts with narrow-band [O m l and Ha images which 
show prominent emission here dBuckalew fc Kobulnickvl 
l200fih . There is also a no n-thermal radio source at this po- 
sition l|Greve et alj|2002h . Furthermore, CO emission is ob- 
served at this location l|Tavlor et all 1 19991) . possibly indi- 
cating this s ource is an embedded, youn g star cluster, as 
suggested bv lWestmoquette et al. (l2007bl ) . In the ir detailed 
study of this area, Westmoquette et al. (|2007bl ) were not 
able to detect any stellar signatures. Thus, our detection 
of [Hei] is the first proof that this source is a young clus- 
ter forming from a blister in an expanding bubble of gas 
i| Westmoquette et alJkOQTbl ). 

Just south of SSC A appears another region of emission. 
It is most prominent in [Hei], but also seen in Pa/3 and 
continuum. This source is a known Hn region JWallerlll99ll) 
and is al so seen in [O m] images (|Buckalew fe Kobulnickvl 
I2006T I. In I Westmoquette et"ai] (|2007bh . there is evidence for 
an unusual clump with a high velocity and large line widths. 
These authors suggest this clump lies along our sight-line 
and is not necessarily connected with the central regions 
of the galaxy. In our equivalent width map, we do not see 
evidence for large line widths here, but then we lack the 
spatial resolution of CMOS. 

One of the most intriguing parts of NGC 1569 sits 
NW of SSC A, where star formation still appears to be 
going on. Optical images show an abundance of star clus- 
ters a nd bright stars here l|Hunter et al. 200d:lAnders" et al.l 
|2004| ) . Many of the clust ers in this region have ages of ~25 
Myr (jAnders et al. | |2004|V The region is dom inated by clus- 
ter 10 as designated by Hunter et al.l ((2000). Narrow- band, 
He 11 observations show emission coincident with the cluster, 
and the authors hypothesize this e mission is due to thre e 
WNL stars or is nebular in origin l|Buckalew et al.ll2000l ). 
Recent, optical luminosity measurements of cluster 10 con- 
firm that massive stars are the bes t explanation for this 
emissio n (|Westmoquette et al.ll2007al ). IWestmoquette et al.l 
l|2007al ) acquired HST images of this cluster and show that 
it is actually two clusters, 10A and 10B. Using photome- 
try, they estimate 10A is between 5-7 Myr, w hile 10B is <5 
Myr. In addition, cluster 6 (as design ated by Hunter et al.l 
(2000)) also has a young age, 4 Myr (| Anders et al.ll2004) . 
These cluster ages indicate a younger po pulation of clusters 
here. Furthermore, iPasquali et al.l (|201ll ) find evidence for a 
higher star formation rate in the NW. 

South of cluster 10, and on the edge of thi s large 
[He i] complex, sits the Hn region identified by IWallerl 
(1991). This nebular region is coincident with [O ill] 



dBuckalew fc Kobulnickvl |2006h and [Siv ] (jTokura et all 
I2006T ) emission. It also sits perched on the edge of a 
giant molecular c loud (CMC) seen in CO observations 
(Taylor et allll999l ). 

Progressing farther south we find the strongest peak in 
both [Hei] and Pa/3. The Pa/3 peak is offset to the north of 
the [Hei] peak by 0.6 arcsec (see Fig. 1). Considering the 
average shifts in the image slices used to create the recon- 
structed image is 0.3±0.8 arcsec, we can not reliably differ- 
entiate the positions of the [He 1] and Pa/3 peaks. 

Interestingly, the [Hei] peak correspo nds with the 
bright est mid-IR, [Siv] peak observed by iTokura et al.l 
(2006). These authors suggest this is an embedded cluster 
consisting of 40 07 stars. Our observations of coincident 
[He 1] emission seems to confirm their finding and is the first 
conclusive evidence that a cluster is forming here. Further- 
more, comparing the locat ion of the [He i] an d [S iv] peaks 
to CO maps presented in jTavlor et al.lll999T ), suggest this 
cluster is buried in the CMC that extends to the southwest. 

This discussion demonstrates how newly formed star 
clusters provoke the continued formati on of star clu s ters in 
the local environment of NGC 1569. lAnders et all (|2004T ) 
mention that star formation is occurring along bubble walls, 
indicating self-propagating star formation, which isdis- 
cussed by several authors. iBuckalew fc Kobulnickvl (|2006T l 
find that the surface density of Hn regions remains con- 
stant with in annuli as the annuli are moved away from 
the center. They suggest this indicates feedback effects are 
confined to the cen t er of the galaxy on scales of 1 pc. 
IWestmoquette et al.l (|2007bl ) use HST observations to show 
the cometary appearance of molecular clouds, indicating 
the effects of stellar winds from the inner clusters. Lastly, 
IPasquali et al.l (120111 ) find that line ratios of [Fell] 1.64/xm / 
Br7 are highest at the edge of gas-filled cavities and in the 
outskirts of the galaxy. This indicates the average between 
past and on-going star formation. In addition, it indicates 
that star formation has been quenched in the interior and 
more recently, triggered in a ring around the cavities ~4 
Myr ago, after the SSCs formed. Our observations of [Hei] 
coincident with CO clouds provide evidence for embedded 
clusters that were triggered to form by cluster winds near 
the CO clouds. 



4 CONCLUSIONS 

Our observations of clumped [Hei] emission to the N of SSC 
A and in the NW provide the first, conclusive evidence of 
newly forming star clusters here. Combining our observa- 
tions, with additional, multi-wavelength data, we discussed 
a scenario of the on-going star formation in this galaxy. 
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